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Abstract. This paper considers edge interval estimation between two regions of a Syn- 
thetic Aperture Image which differ in texture. Different point estimation strategies under 
multiplicative noise are discussed in the literature. It is important to assess the quality of 
such point estimates and to also perform inference under a given confidence level. This 
can be achieved through interval parameter estimation. To that end, we propose bootstrap- 
based edge confidence interval. The relative merits of the different inference strategies are 
compared using Monte Carlo simulation. We also analyze a real dataset. 



1. INTRODUCTION 

Imagery obtained through the use of coherent illumination suffers from a noise known 
as speckle. This is the case, for instance, of synthetic aperture radar (SAR) images. The 
noise is usually modeled in a multiplicative fashion and different inference strategies are 
used in order to extract useful information from data subject to such speckle noise. Much 
attention has been devoted to statistical inference using SAR data because the use of SAR 
sensors has become widespread. Such sensors are particularly useful because they do not 
require external sources of illumination and their wavelength is not affected by weather 
conditions. There are three possible types of texture for each region of a SAR image, namely: 
homogeneous (i.e., surfaces with little texture such as lakes, deforestation areas, crop fields, 
deserts and snow covered areas), heterogeneous (i.e., surfaces with more texture than ho- 
mogeneous areas, but with not as much texture as forests at regular relief) and extremely 
heterogeneous (i.e., areas that have very intense texture, such as urban areas). 

Several approaches for detecting an edge between neighboring regions with different tex- 
tures have been proposed in the literature. Gambini et al. |1| proposed a parametric point 
estimator (i.e., detector) for the edge location. It has been shown that their estimator typ- 
ically outperforms alternative detectors I21|31ID in a number of different settings |5|. Non- 
parametric detectors have also been proposed. Bovik et al. |6| used nonparametric meth- 
ods to detect edges under additive Gaussian noise. Beauchemin et al. |7| used an alterna- 
tive approach based on the Wilcoxon-Mann-Whitney statistic to detect changes in adjacent 
sets of pixels. Lim and Jang 18j used two sample tests to detect edges in images subject 
to noise. Giron et al. [?] introduced detectors based on the Mann-Whitney, squared ranks 
to variances, Kruskal-Wallis and empirical distribution nonparametric test statistics. The 
authors used Monte Carlo simulation to compare the small sample performances of their 
detectors to that of Gambini's estimator. Their numerical evidence shows that the Kruskal- 
Wallis and Gambini's detectors behave similarly, the former however being considerably 
less costly from a computational viewpoint. Indeed, their numerical results show that edge 
detection using the Kruskal-Wallis detector is approximately one thousand times faster than 
that based on Gambini's maximum likelihood detector. 
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Figure 1 . Rectangular window with ortiiogonal edge to line detection. 

All edge detection strategies that have been proposed in the literature yield edge point 
estimates. Our main goal in this paper is to consider interval edge estimation. We intro- 
duce different strategies that can be used to produce a confidence interval for an edge that 
separates two regions of a SAR image v«th different textures. The interval estimators are ob- 
tained by setting the desired confidence level and their lengths are indicative of how certain 
one can be that the edge has been located. 

More specifically, we wish to detect the position of the edge between two regions along 
a line segment (detection line) . To that end we use information on the pixels that lie in a 
neighborhood of a rectangular shaped detection line (window). Figure [T] shows a simulated 
SAR image that contains two distinct regions and also a detection window divided by the 
detection line. 

As pointed by Gambini et al. (H, their detector can deliver poor estimates when the de- 
tection line is not orthogonal to the boundary, as it is the case in Figure[T] In contrast, edge 
interval estimation does not rely on the orthogonality between the detection line and the 
boundary. This reinforces the importance of interval edge detection. 

We chose to construct confidence intervals based on the Kruskal-Wallis point estimator 
because the available numerical evidence shows that it performs as well as the maximum 
likelihood detector (Gambini's detector) at a much lower computational cost. We also con- 
sider different bootstrap interval estimators: normal approximation, basic bootstrap, per- 
centile and two variants percentile- f method that we propose. For details on bootstrap in- 
terval estimation, the reader is referred to Efron and Tibishirani 1 9 1 . Our numerical results 
show that the bootstrap percentile estimator is the overall best performer in the sense that 
it typically has the best coverage amongst all estimators. 

The remainder of the article unfolds as follows. Section|2]introduces the statistical model 
we shall use to model SAR data. Section [s] presents the Gambini |1, 5| and Kruskal-Wallis 
detectors. Sectionjijconsiders edge bootstrap interval estimation. Sectionjsjpresents Monte 
Carlo evidence. Sectiorjejcontains an application that uses real (not simulated) data. Finally, 
Section[7]offers some concluding remarks. 

2. The SAR image model 

The multiplicative statistical model is widely used in SAR image data analyses. According 
to this model, the data are described by a random variable Z which can be viewed as the 
product of the independent random variables X and Y, where X models the properties of 
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Figure 2. ^1 densities for L e {1, 3} and a e {-2, -5,-15} (dots, daslies, and solid). 



the imaged area (backscatter) and Y models the multiplicative noise (speckle) introduced 
by the use of coherent illumination which degrades the image quality. 

Images represented in the intensity format can be described by the distribution de- 
noted by Z ~ 'S^. The parameters that index such a distribution are: (i) the number of looks 
L> I, which is a measure of the signal-to-noise ratio, (ii) the scale parameter 7 > 0, which 
is related to the relative strength between the incoming and reflected signals, and (iii) the 
roughness parameter a < 0, which relates to the land type texture. The larger the value of 
a, the more heterogeneous the area: when a < -10 the area is very homogeneous (e.g., pas- 
tures), when -10 < a < -4 the area is heterogeneous (e.g., forests) and when -4 < a < the 
area is extremely heterogeneous (e.g., urban areas). In what follows we shall assume that 
the number of looks is known. 

According to the model proposed by Frery et al. , the speckle noise is gamma dis- 
tributed, denoted by F ~ T{L, L), its density being 

tL 



fYiy) ■ 



T[L) 



^exp(-Ly), I> l,y>0. 



(1) 



The backscatter can be modeled using the inverse gamma distribution, X~T ^ (a, 7), whose 
density function is 

1 



/x(x) = 



-a-l 



7«r(-a) 



exp 



I- }• 



a,7,x > 0. 



Thus, Z = XF is distributed as Z ~ '^^{a,7,L), whose density function is 



fz{z) 



L^T{L- 



a) 



,L-1 



L> l,-a,7,z> 0. 



(2) 



(3) 



7«r(I)r{-a) ij + Lz)^-" 

Figure[2]shows the density for L e {1, 3} and a e {-2, -5, - 15}. 

It is noteworthy that random number generation can be easily performed. It suffices 
to use the fact if ~ T~^(_k,d) then l/W ~ T(_k,l/6). Thus, in order to generate a pseudo- 
random number from we only need to generate a pseudo-random number from T{L, L) 
and then divide it by a pseudo-random number obtained from r(-a, I/7). 
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Using l[3j, it can be easily shown that the rth noncentral moment is 



ElZ'-] = [^] 



JY ri-a-r]T[L+r) 
Ti-a)T[L) ' 



(4) 



if -a > r, and oo otherwise. Notice that the value of 7 that corresponds to ElZ] = 1 is 

n-a)TiL)L 

T= • (5) 

' T{-a-l)T{L+l) 

When a < -I, parameter estimation can be performed using the method of moments. 
Let {Zi,...,Z„) be a vector of independent and identically distributed random variables with 
common distribution (a.y.L). The rth sample moment is given by 



mr = -f^z'.. (6) 



Using jijl we can obtain expressions for the moments of order | and 1 as a function of a 
and 7. Replacing the population moments by their sample counterparts and replacing the 
parameters a and 7 by their respective estimators d and 7, we obtain the following system 
of equations: 



mi 



mi 



'7] r(-d-l)r(L+l) 
.ij T{-a)T[L) 

(7) 

'yY r(-d-i)r(L + i) 



\L) n-d]T{L) 



It then follows that 



miT{-a)T[L)L 1 



r(-d-i)r(L+i) YH-d-lwHL+l)' 

which can be numerically solved. After obtaining the moments estimate of a we can then 
compute 7. 

In what follows, we consider that there is an edge when it is surrounded by two regions 
with textures. In this case, one of the regions is represented by ^^ia^,Ye,E) whereas the 
other region is represented by ^^{ar,Yr>L). 

3. Edge detection in SAR imagery 

3.1. The Gambini parametric detector. The Gambini algorithm is based on the fact that if a 
point belongs to the edge, then a sample taken from its neighborhood is expected to exhibit 
a noticeable change in the parameter values which are used to describe the pixel distribu- 
tion on either side of the edge. Let 5 be a line segment (detection line) such that a point 
that belongs to 5 also belongs to the edge that separates two regions. The interest lies in de- 
termining which point detection line is a transition boundary. The segment 5 = (zi, . . . , Zjv) 
is a strip of pixels obtained by discretization of a straight line on the image. Assume that 
there is a nonnuU intersection between the edge and the detection line and let zj be a point 

such that [zi,...,Zj) comes from the '^^(a^, 7^, L) distribution and (zy+i Zjv) comes from 

the '^^(ar,7r,i) distribution. Then, the value of j indicates the edge location. To find the 
transition point at 5, consider the likelihood function given by 

L = L(_a(,Ye,ar,Yr)=WPr{Zk;a(,Y()>^ H Pr^Zk;ar,Yr)- (8) 
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The log-likelihood function 

/ N 

= lnL= J^lnf^o[Zk;a£,re)+ ^ In f^^o {Zk) a r.jr), 

k=l k=j+l 

where j e {l,...,Ar- 1}, is maximized when the index j is the edge. Using Equation l[3|, it 
follows that 



^= y In-r 

k=\ f"/T{L)T{-af)ife + Lzk)^-^f 

+ y In—. . (9) 

k=j+i f"'T{L)T{-ar){fr + Lzk)^-"r 

Thus, the estimated transition point at 5 is given by 

jGE = argmaxi!?. (10) 

The Gambini's estimator Jge requires the estimation of (a^,7^) and [ar,Jr) for all possi- 
ble values of j e {l,...,N -I}. Such estimations are carried out using rectangular windows 
in which the detection line coincides with the major axis of the window. Note that for ex- 
treme values of j estimation is performed using a small sample, and as a consequence the 
resulting estimate can be poor. We also note that the moments estimator can be quite bi- 
ased. 

3.2. Nonparametric edge detection: Edge detection using nonparametric methods as al- 
ternative to the Gambini's estimator was proposed by Giron et al. [?] . In what follows we 
shall use their best performing estimator, namely: Kruskal-Wallis. 

Consider k independent random samples (with possibly different sizes) from k identical 
populations or from k populations such that at least one of them tends to produce obser- 
vations v«th larger values, the ith sample containing n observations. Let N = Y.\^.^ni be the 
total number of observations. The Kruskal-Wallis statistic was developed to test the null hy- 
pothesis that the fc-samples come from the same population. If the null hypothesis is false, 
then at least one of the k samples will tend to produce observations with larger values. In 
order to compute the Kruskal-Wallis statistic we have to consider N observations and assign 
ranks to them. Let R{Zij] be the rank of Zij, the jth observation of the ith sample, and let 
Ri = Rizi j) be the sum of ranks to the ith sample. The test statistic is given by 



I (^Rj jV(iv+i) M 



where 



1 



N-l 



LLRi^ij) T- 

\i=lj=l * 



If the number of ties is sufficiently small (or null), it can be shown that the statistic T can 
be expressed as 

12 ^ 

Tkw= ^ y — -3(jV+l). (11) 

N{N+ 1) rii 
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Edge estimation using Tkw is straightforward. Consider a strip of pixels s = (zi zn) 

and assume that the jth element is the edge. The edge thus splits the data into two sam- 
ples, (zi, . . . , Zj] and (zj+i,. . . , zn), whose sample sizes are «i = j and «2 = N-j, respectively, 
where 7 £ {1, . . . , iV- 1}. We conclude that the two samples are come from different popula- 
tions when Tjcw is large, the corresponding edge estimator being 

jKw = argmax Tkw- (12) 

j 

This edge estimator is computationally less costly than Gambini's parametric estimator. 



4. Bootstrap-based edge confidence intervals 

It is important for practitioners to be able to compute not only edge point estimates but 
also edge interval estimates, that is, confidence intervals for SAR image edges. We note 
that the distribution of KW point estimator is unknown and hence the usual corresponding 
confidence intervals cannot be obtained. An alternative is to use bootstrap resampling to 
construct interval estimates. 

In what follows we shall consider bootstrap edge interval estimation. In particular, the 
following bootstrap interval estimators shall be considered: normal approximation (NAP), 
bootstrap basic method (BBM), percentile (PERC), and Studentized (or percentile- boot- 
strap- t) (ST). Given the high computational burden of the percentile- 1 method, we propose 
two new variants of such a confidence interval that are much less computer-intensive. 

Consider again a strip of pixels s = (_Zi,...,zjs[), with an edge between two regions lo- 
cated at j e {l,...,N}. Let j]cu, be the KW edge estimate. Using the estimate jf^u/ we can 
obtain pseudo-samples s* = (zJ' , . . . , zj^) , where (z^ , . . . , zt ) is obtained by randomly sam- 

./ k w 

pling with replacement from (zi , . . . , z ^ ) . Similarly, (z? , . . . , zt,) is obtained by randomy 

Jkw Jkw~^^ 

sampling with replacement from (zj^^^j,...,ZAr). Let j* be the edge estimate obtained by 
resampling a strip of pixels s*. After executing this scheme B times, we obtain B edge es- 
timates (jj ,...,7g). This subsampling process is known as nonparametric bootstrap. Let 
j*iy ■ ■ y 7(*B) the ordered values of jy,..., j^. The empirical distribution of Jkw is given by 

#{/* < /} 

P*bU') = ^^ be{l,...,B}. (13) 

It can be used for constructing confidence intervals for the true edge location. 

Assuming that Jkw is a consistent estimator of the edge and given that the true distribu- 
tion of Jkw is unknown, a simple approach to construct confidence intervals is to assume 
that the distribution of Jkw - j is N{0, v), where v is the variance of Jkw- One can estimate 
V by the sample variance of j*, i.e.. 



v = — 
B 



1 ^ f^* 1-^2 



where /* = {XI B^Y^^xilby LetO< a< l/2andzi-ai2 = ^ kl- a/2), i.e., zi-a/z is the l-a/2 
standard normal upper quantile. A 100(1 - a)% confidence interval for the edge location is 

Jal2'Jl-al2^NAP=ljkw-V Zi-al2,]kw+V Z\-al2\, (14) 
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We can improve the interval in (14) by bias correcting the edge point estimator. The bias 
of jkw obtained from strip 5 can be estimated as 

^* = ^* Ukw) = ]* - hw 

Thus, after bias correction, the bootstrap confidence interval with coverage 100(1 - a)%is 

[j*al2'Jl-al2^NAP = Ijkw-^* " V^'^ Zi-all, Jkw - ^* + V^'^Zi.al2] ■ (15) 

It is also possible to construct an edge confidence interval based on the bootstrap ap- 
proximation to the distribution of Jkw - j- Here, the quantiles are obtained by using the 
ordered values of j^j^ - j^^jj. The 100(1 - a)% basic bootstrap method (BBM) confidence 
interval is 

[jal2'j\-al2^BBM = [jkw - ^i*B(l-al2)) ~ jkw): jkw - U*B(a/2)) ~ jkw)] ■ (16) 

An alternative bootstrap confidence interval is the Percentile interval. The 100(1 - a)% 
percentile interval is given by 

[jll2'j\-al2^PERC = U*B(al2))'j*B{l-al2))] ■ ^^'^^ 

It is noteworthy that this interval is may be asymmetric and only contains proper values of 
the edge. 

An alternative approach is to construct a bootstrap studentized confidence interval: the 
percentile- 1 or bootstrap-? interval. To that end, in each bootstrap replication we compute 
the quantity 

where vl^'^ is the standard error j* obtained in the bth bootstrap replication. The ath sam- 
ple quantile, , is obtained as follows: 

#{Z; < \a)} 

= a, 

B 

b e {I,... B}. The 100(1 - a)% studentized bootstrap confidence interval is then given by 

[]al2>]l-al2^ST = ijkw " V^'^kB{\-al2)),}kw " V^'^hBiall])] ■ (19) 

The sample variance v*^ can be computed by using a sub-bootstrap (i.e., a second level 
bootstrap) with B' replications, where B' = 50 is the minimum number of resamples that 
yields reasonable estimates. Efron and Tibshirani |9J recommend using B' = 200. However, 
they also recommend using B = 1000. This implies that estimation of v* using a second 
level bootstrap adds considerable computational burden, even when one uses B' = 50. 

In what follows we shall propose two alternative approximate bootstrap- f interval esti- 
mators, denoted STl and ST2. By using them one can compute v* without having to resort 
to a second level bootstrap. 

The idea behind out STl bootstrap method is to estimate i^* using a subset of size B' from 
the replications 7\* , . . • , jg imposing the restriction that the sample variance is positive. Let 
j*' = {Ji > • • • ' jg'} be the randomly selected subset of B bootstrap estimates. Our algorithm 
can be outlined as follows: 

Algorithm 1: STl Method 

(1) If 7^ = jkw< then Z* -0 (see Equation (18)) and the algorithm ends. 

(2) If the variance ]{,■■■, jg is not equal to zero: 
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(a) Obtain j* - {j^ j*, ] by drawing with replacement from 

(b) Compute 



1 ^' 



( . I B' 



^ ^ b'=l \ ^ b"=l 



\ 

(c) If - 0, go back to step 2{a). Otherwise, the algorithm ends. This step is to be exe- 
cuted at most B" > times. 

(d) If the limiting number of repetitions for step 2(c) {B" ) is reached, set y^ equal to the 
sample variance of yj* , . . . , jg ■ 

(3) If the sample variance of jj* , . . . , 7^ is equal to zero: 

(a) Obtain 7 * by estimating the edge using a new bootstrap replication of pixels strip s. 
Repeat this step until j* ^ . This step is to be executed at most B" > times. 

(b) If/:' = do = + 1. 

(c) Set equalto the sample variance of {7j,...,7*,_j, 7* }. 

It is noteworthy that Algorithm 1 must be executed for each be{l,...,B}. This algorithm 
is considerably less computationally intensive than the standard bootstrap-? algorithm in 
which a second level bootstrapping scheme is used for variance estimation. 

We shall now introduce a second approximation to the bootstrap- 1 method. 

Algorithm 2: ST2 Method 

(1) Obtain estimates of the edge based on Bx new resamples from the pixels strip s, thus ob- 
taining {j*^ ^ J'^^J . 

(2) If the sample variance of the values 7 * is not equal to zero: 

(a) Obtain 7* - [jf ,. . . , 7'^,} [B' < B^) by drawing with replacement from 7* . 

(b) Set equal to the sample variance of 7* . 

(c) If - 0, go back to step 2(a). Otherwise, the algorithm ends. This step is to be exe- 
cuted at most B" times (0 < B" < Bx). 

(d) If the limiting number of repetitions set for step 2(c) [B") is reached, set i)^ equal to the 
sample variance of 7 * . 

(3) If the sample variance of the 7 * values is equal to zero: 

(a) Obtain 7*' by estimating the edge using a new bootstrap replication of pixels strips s. 
Repeat this step until that 7*' ^ Jb+v This step is not to be executed more than B" 
times. 

(b) ifi;' = iB+i'do/:' = /B+i + i- 

(c) Estimate 1)^ by sample variance of {7^+p . ■ . , ig+B'-vi* 

Notice that STl is less computationally expensive that the ST2 and that they are both 
less computationally costly than the standard bootstrap-?, in which a second level (inner) 
bootstrap is carried out for variance estimation. 

5. Numerical RESULTS AND DISCUSSION 

In what follows we shall report the results of several Monte Carlo simulations which were 
performed to assess the finite sample merits of different edge interval estimates in SAR im- 
ages. All simulations were carried out using Ox matrix programming language (1 11 1). The 
hardware used was an Intel(R) Core(TM)2 Quad CPU Q6600 2.40GHz computer running on 
Ubuntu Linux. Graphics were produced using R programming environment 1 12 1 . 

In all simulations, L = 1. Recall that this is the most challenging situation. All images 
we generated are rectangular and contain 20 x 100 pixels {N = 100) with an edge located 
at j = 50. Thus, each image contains two distinct regions of equal size. In the region to 
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the left of the edge a = a(, and otherwise a = Ur, with a(,ar e {-2, -3,..., -15}. For each 
region, the scale parameter y was chosen according to |5}. The most challenging situation 
as far as edge detection is concerned takes place when |ar - a^l = 1; the higher the absolute 
difference between and a^, the easier it is to locate the edge. 

The number of Monte Carlo replications is R = 5000 and the number of bootstrap repli- 
cations is B = 1000. For each configuration, we constructed 5000 confidence intervals using 
the normal approximation approach, the basic bootstrap method, the percentile method, 
and also the STl and ST2 methods. The nominal coverage is 1 - a = 0.95 (95%). The results 
for the STl and ST2 methods were obtained using B' = 50, and B" = 200, and for ST2 we 
used Bjc = 200. 

We computed the exact coverage (C) of each interval estimator corresponding to the 95% 
nominal level, which was done for each configuration {af,ar). The best performing esti- 
mator is the one whose exact coverage is closest to the nominal coverage. Figure|3]plots the 
distance D = \C- 0.95| x 100% against a,- (the roughness parameter value to the right of the 
edge). The six panels in the figure correspond to a£ e {-2, -5, -8,-11,-13,-15}. Situations 
in which ag = ar (i.e., there is no edge) are not considered. 

Visual inspection of Figure |3] reveals that all interval estimators perform equally well 
when ag = -2 and < -4; their coverages is 100%. We also notice that as - a^l in- 
creases the differences between the different coverages tend decrease, i.e., the different in- 
tervals display similar coverages. That happens because the distribution of Jkw tends to be 
more concentrated around of the true value j = 50 as - a^l increases, i.e., the estimator 
Jkw is more accurate when the two regions have very distinct textures. Such a tendency is 
not, however, uniform. Also, when the regions on both sides of the edge are very heteroge- 
neous (for example, a£ = -2, and = -3, or ae = -5 and a£ = -4), the different intervals 
behave similarly, but when the two regions of are very homogeneous (for example a£ = -13, 
and ar = -12), the percentile estimator outperforms all other interval estimators. Figure|3] 
shows that the percentile estimator is the best performing interval estimator in most situa- 
tions. It is also noteworthy that the percentile estimator is the best performing estimator in 
the most challenging situations, i.e., whenever the two regions are very homogeneous and 
have similar textures. 

In Figure [ij a symbol indicates which is the best performing interval estimator for each 
configuration [a^, a^); the quantity bellow the symbol is 

^ ~ ^percentile ~ ^best method' 

which measures the method performance relative to the percentile method. The symbols 
used in the figure are the same as the ones used in Figure |3] except for the symbol used 
to represent BBM, which is now is>. The symbol M is used whenever all methods produced 
the same coverage. We observe that in configurations in which the two absolute values of 
roughness parameter are large and similar (the most challenging situations), the percentile 
estimator is the best performer, closely followed by BBM. 

The two studentized estimators, STl and ST2, outperformed the competition whenever 
\ae - ar \ < 5. In these situations, the values of A are between 1% and 2%. Since these two 
estimators behave similarly, we recommend the use of STl, which is less computationally 
intensive. It should be noted that in such situations the percentile estimator is quite com- 
petitive vnth the studentized estimators. 
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Figure 3. Distances between empirical and nominal coverages. 



The normal approximation interval estimator becomes the best performer when \ae-ar\ 
is large, being closely followed by the percentile estimator. All methods behave similarly 
when the absolute difference between a/ and becomes very large; see Figure|3j 
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Figure 4. Best method and difference performance A between the per- 
centile and the best performer. 



When taken together, our numerical results show that the percentile estimator is either 
the best performing estimator or quite competitive with the best performer. Therefore, we 
conclude that the percentile edge estimator is to be preferred. 

We have also carried out simulations in which there is no edge, i.e., a£ = ar- Here, we 
shall focus not on coverages but on the average lengths of the different confidence intervals. 
The estimated edge locations are expected to be randomly and uniformly distributed along 
the detection line since there is no edge. We thus expect the average interval length to be 
approximately equal to the detection line length. Figure [s] displays the the average interval 
lengths. We notice that all methods produce similar average interval lengths, regardless of 
the configuration (a^, a, ), as expected. The percentile and BBM interval estimators have the 
same average length: approximately 93 pixels. The STl and ST2 estimators displayed similar 
average lengths (approximately 95 pixels), the former always being slightly narrower. The 
average lengths of the normal approximation interval always exceeded the detection line 
length (100 pixels) . Overall, we notice that all methods yield very vnde intervals when there 
is no edge to be located. A wide confidence interval can then be taken as evidence that there 
is no edge. 
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Figure 5. Average length for situations where there is not edge. 



Figure [e] displays the average interval lengths. Notice that when \ a^-ar \ is large, the av- 
erage lengths are very small. On the other hand, when | - a ^ I is small, the average lengths 
tend to be large. Normal approximation yields the smallest average lengths whenever the 
average lengths are not all the same. The next best performers are the percentile and BBM 
intervals. The largest average lengths correspond to the studentized interval estimators, STl 
and ST2. The largest difference in average lengths (7 pixels) is obtained when we consider 
the NAP and percentile intervals. The average lengths of the studentized and percentile 
intervals are very close. 

The execution times of the different interval edge estimators for each Monte Carlo repli- 
cation are given in Table[T] Since the computational cost does not depend on the parameter 
configuration, we use (a^, a^) = (-2, -3) in order to measure the different execution times. 
We report both the total time (in seconds) it takes to compute the interval estimator and that 
relative to the bootstrap- 1 estimator (%). The bootstrap- f estimator was computed using 50 
replications in the inner bootstrap, which is carried out for variance estimation. The results 
show that the Bootstrap- 1 is approximately 50 times more computationally intensive than 
all competing methods. The normal approximation, BBM, and percentile methods have 
nearly the same computational cost. STl and ST2 are slightly more costly, STl being less 
computationally intensive than ST2. It should be noted that the percentile method is not 
only the best performer as far as coverage is concerned, but it is also the least computation- 
ally intensive method. 

6. An APPLICATION 

In this section we apply the different edge interval estimators to a real image, i.e., to ob- 
served (not simulated) data. The SAR image is displayed in Figure [Tj The rectangles in 
the image are the regions chosen for edge detection. Rectangles (regions) 1,2 and 4 have 
210 X 100 pixels whereas rectangle (region) 3 has 100 x 210 pixels. In the edge detection pro- 
cess, we divided rectangle 3 in ten horizontal windows, each of 21 x 100 pixels, the detection 
line being horizontal for each window and dividing the vnndow in equal parts of 10 x 100 
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Figure 6. Average interval lengths for ag £ {-3, -8,-13} (from top to bottom). 
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Table 1 . Computational cost and computational cost relative to bootstrap- 1. 



Method Execution time (seconds) Relative to bootstrap- 1 (%) 



Bootstrap-f 81.422 100.00 

NAP 1.594 1.96 

BBM 1.578 1.94 

STl 1.630 2.00 

ST2 2.191 2.69 

PERC 1.565 1.92 




Figure 7. E-SAR image with 1100 x 2400 pixels, and delimited regions. 



pixels. Similarly, each of the other rectangles was divided in ten windows of 100 x 21 pixels, 
the detection line being vertical and dividing the window in equal parts of 100 x 10 pixels. 
Therefore, for each rectangle we compute ten interval edge estimates. In the figures that 
follow, we use green lines to connect the lower interval limits and also the upper interval 
limits; red lines are used to connect the different edge point estimates. Each detection win- 
dow is delimited by yellow lines, the windows being numbered from left to right or from top 
to bottom (for example, the first window in rectangle 1 corresponds to the window located 
in the extreme left). 

Figure [s] contains results relative to edge detection in rectangle 1. We observe that the 
confidence intervals are quite narrow (the smallest interval length equals 3 pixels), except 
for the first two detection windows. In the second window there are less pixels in light re- 
gions than in dark regions, in which there is more noise, hence the wider confidence in- 
terval. This happens because edge detection becomes quite challenging under complex, 
rich textures. We also note that the estimated edge is typically located in areas in which the 
differences in texture are largest; see, e.g., the first window. 

A challenging situation for edge detection can be seen in rectangle{region) 2 of Figure 
[yj There we notice two fairly different textures, but it is not easy to identify where exactly 
lies the frontier. The results of edge detection for this region are presents in Figure |9] Note 
that, in all detection windows, the confidence intervals are fairly wide; the smallest interval 
length is 5 pixels (last detection vnndow) and the widest interval covers 26 pixels (next to 



INTERVAL EDGE ESTIMATION 



15 




Figure 8. Edge interval detection in region 1. 




Figure 9 . Edge interval detection in region 2. 



last detection window). In challenging situations such as this interval estimation becomes 
quite useful since it signals that point edge detection may not be reliable. 

Edge detection in rectangle (region) 3 of Figure [t] is even more challenging that in rec- 
tangle (region) 2 because the edge divides the image in two regions with textures that are 
visually very diffuse. The detection results are displayed in Figure 10 We note that the in- 
terval lengths are nearly the same in all detection windows. None of interval lengths exceeds 
1 1 pixels. We conclude that, despite the aforementioned difficulty, the KW detector yields 
accurate edge estimates. 

A limitation of the point edge detection methods we consider is that an edge will be al- 
ways located even when there is none. This happens because T^w in flT) always achieves 
a maximum value. Interval estimation can, however, be used to assess whether an edge in- 
deed exists. Such a situation is explored in rectangle (region) 4 of Figure [Tj A confidence 
interval that covers 94 pixels is obtained, thus suggesting that there is no edge. In the last 
window detection, we obtain a fairly wide interval (41 pixels) which can be taken as evi- 
dence that the point edge estimate is not accurate. A similar situation takes place at the first 
detection window. Here, however, the dark regions are most likely not the result of noise but 
reflect the terrain characteristics. 
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Figure 10. Edge interval detection in region 3. 




Figure 1 1 . Edge interval detection in region 4. 



7. Concluding remarks 

We addressed the issue of edge detection in SAR images. Several bootstrap-based interval 
estimates were described and numerically evaluated. Overall, the percentile confidence in- 
terval proved to be most reliable, especially in the challenging situation in which the regions 
on both sides of the edge have similar textures. The percentile interval typically displayed 
the best coverage but it is typically wider than alternative intervals. There is thus a trade-off 
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between coverage and length. The normal approximation method delivered intervals with 
the smallest lengths. 

We proposed two variants of the bootstrap-? method for edge interval estimation, de- 
noted STl and ST2. They are less computationally costly than the standard Bootstrap-? 
method. They tend to work equally well, STl being more computationally efficient. 

We have considered situations in which edge detection is carried out in an image region 
in which there is no edge. Point estimation will always locate an edge, even when there is 
none. In such situations, the resulting interval estimates tend to be fairly wide, thus sig- 
naling that the detected edge is not to be trusted. Very wide intervals can be taken as an 
indication that most likely there is no edge in that region of the SAR image. 

Finally, real (not simulated) data were analyzed. We performed edge detection in several 
regions of a SAR image. 
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